Topologically non-trivial superfluid phases and Majorana fermions from Kohn-Luttinger effect 
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Spin-triplet p-wave superfluids of single-species (or spin-polarized) fermionic atoms is a topological super- 
fluid. In 2D such a superfluid supports zero-energy topological Majorana fermion excitations in order parameter 
defects such as vortices and sample edges. In 3D these superfluids support topologically protected Dirac points 
in the bulk spectrum and flat surface Majorana arcs. Despite the promise, creating a spin-triplet p-wave single- 
species fermionic superfluid either from direct p-wave Feshbach resonance or indirectly from the recently pro- 
posed artificial spin-orbit and Zeeman couplings can be technologically challenging. Here we show that such 
topological superfluids can be far more simply created by using the Kohn-Luttinger effect applied to two species 
of spin-polarized fermions with a density imbalance. We discuss how the topological Majorana excitations of 
the resulting superfluids can be identified using recently-developed experimental techniques. 



Introduction: Spin-triplet p-wave superfluids of spin- 
polarized fermionic atoms offer an excellent test-bed for ex- 
ploring many intriguing quantum phenomena |ll]-{3 , 47]. Such 
systems are generically topological in all spatial dimensions, 
exhibiting, in 3D, topologically protected Dirac points in the 
bulk spectrum 1 1 , 5, 6] and surface Majorana arcs |7], while in 
£' = 2,1 they exhibit zero-energy Majorana fermions in order 
parameter defects such as vortices and sample edges 12. Si- In 
recent years, the physics of the 2D chiral p-wave (px+ipy) su- 
perfluids has attracted much attention yj because of its non- 
trivial statistical properties ll47ll and potential application in 
topological quantum computation jsilsl-fioll. The chiral super- 
fluids can also act as a testbed for studying quantum phenom- 
ena such the quantum non-locality and the violation of Bell's 
inequality [10], which are often masked by the many -body 
effects in a macroscopic system. Schemes for observing any- 
onic statistics and implementing topological quantum compu- 
tation using vortices in these systems have also been proposed 
llO.l . A potential advantage of the spin-polarized cold atom 
Px + ipy chiral superfluid over solid state chiral-p wave sys- 
tem such as strontium ruthenate 111 llll2[l is that, the cold atom 
system being spin-polarized, even ordinary vortex excitations 
of the superfluid can have trapped Majorana modes while in 
the solid state system the relevant topological excitations are 
the exotic half-quantum vortices L13ti42J. 

Despite the promise, creating a spin-triplet p-wave super- 
fluid of single-spin (or single-species) fermionic cold atoms 
has turned out to be a challenging task. With the observa- 
tions of p-wave Feshbach resonances in spin-polarized '"'K 
and ^Li atoms in optical traps liTjljlTll . a Px+ipy superfluid of 
fermionic atoms was thought to be soon realizable using Fes- 
hbach resonance [18-20]. However, because of the short life- 
times of the p-wave pairs and molecules fllTll . realizing a spin- 
polarized p-wave superfluid in this route now seems difficult. 
An alternative scheme, involving an artificially created spin- 
orbit coupling and a Zeeman field along with a regular s-wave 



Feshbach resonance, has also been proposed 0211 12211 . Here 



even though the fermions are not spin-polarized, a judicious 
combination of the spin-orbit coupling, Zeeman field, and the 
fermion density can still be obtained to enforce a single Fermi 



surface, thus allowing the formation of a topological super- 
fluid with Majorana fermions. This strategy, which has now 
been theoretically applied also to condensed matter systems 
lE3h28il . still requires artificial laser-induced spin-orbit and 
Zeeman couplings for fermionic atoms which are presently 
beyond experimental reach. Considering the importance of a 
p-wave superfluid of the same species of atoms both for fun- 
damental as well as technological reasons, it's important to 
investigate other, possibly simpler, methods of realizing such 
systems in optical traps without requiring either p-wave Fesh- 
bach resonance or artificial spin-orbit and Zeeman couplings. 
In this paper we achieve this goal by utilizing a completely 
different (and more practical in cold fermions) scheme for 
realizing a triplet p-wave superfluid of the same species of 
atoms by using the celebrated Kohn-Luttinger effect |29-'3ll 
modified to apply to the case of two types of (individually 
spin-polarized) fermionic atoms. The essence of the Kohn- 
Luttinger effect is superfluid pairing in a short-ranged repul- 
sive fermionic system mediated by another fermionic back- 
ground which renormalizes the bare repulsive two-particle in- 
teraction vertex. By using two distinct species of individually 
spin-polarized fermions with a density imbalance we show 
here that it is possible to create single species p-wave topo- 
logical superflids which can be individually probed to uncover 
the topological Majorana excitations. 

Once a single-species spin-triplet superfluid is realized in 
experiments, one natural and important question is how to 
experimentally observe the topological properties, in partic- 
ular, the topological excitations known as Majorana fermions. 
Majorana fermions are zero-energy excitations of a 2D chi- 
ral {px + ipy)-wave or a ID p^^-wave superfluid which are 
trapped by order parameter defects such as vortices (in 2D) 
and sample edges (in ID). Majorana fermions, first proposed 
by E. Majorana to describe neutrinos |'32ll, are quantum par- 
ticles which are their own anti-particles [33] (i.e., the second 
quantized operators 7 satisfy 7^ — 7), unlike Dirac fermions 
where particles and anti-particles (holes) cannot be identified. 
A cold-atom demonstration of Majorana fermions in itself is 
a remarkable experimental goal, but also they have attracted a 
lot of recent attention because of their potential use in topo- 
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logical quantum computation (TQC) yl l8|] . In the condensed 
matter systems Majorana fermions have been proposed to be 
identifiable in local STM tunneling experiments mM and 
in fractional Josephson effect 1271 L 
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However, such 

methods cannot be directly used to detect Majorana fermions 
in cold atom neutral superfluids, because superfluids, in con- 
trast to superconductors, do not directly couple to an elec- 
tric potential difference. Nevertheless, we show how Majo- 
rana fermions in the 2D cold atom spin-triplet superfluids pro- 
duced in our scheme can be identified using radio-frequency 
(RF) spectroscopy experiments. In 2D, according to our the- 
ory, only the majority component (say K) of the two species 
of fermions becomes superfluid in the experimentally achiev- 
able temperatures, while the minority component (say Li) re- 
mains uncondensed. This, however, does not pose any prob- 
lem in imaging the MFs in the majority component since in 
our RF spectroscopy based imaging the individual species of 
atoms can be probed independent of the other species. The 
absence of inter-conversion between the K and Li atoms en- 
sures that the MFs in the majority component remains robust 
even if the minority component remains gapless and uncon- 
densed. Note that a similar robustness of MFs to the presence 
of gapless electrons does not exist in condensed matter sys- 
tems, and, therefore, cold atom fermionic mixtures are partic- 
ularly suitable for realizing MFs which can be robustly probed 
using RF spectroscopy. In the 3D spin-triplet p-wave superflu- 
ids produced in our scheme the topological bulk Dirac modes 
have a structure similar to that in the recently proposed spin- 
orbit coupled 3D superfluid in Ref. [6J. We discuss how these 
non-trivial excitations can be probed in the cold atom ana- 
log of the angle-resolved photoemission experiments. Taken 
together, the method to create a single-species spin-triplet p- 
wave superfluid using the Kohn-Luttinger effect, coupled with 
the methods to observe the non-trivial topological excitations, 
gives a complete description of a promising new way to create 
and analyze a topologically non-trivial superfluid in both two 
and three dimensions. 

Kohn-Luttinger effect and superfluid critical temperatures: 
The essence of the Kohn-Luttinger effect is the possibility of 
a superfluid pairing in fermionic systems with short-range re- 
pulsion that does not support the conventional BCS-type pair- 
ing ll29ll . The presence of a fermionic background renormal- 
izes the two-particle vertex, and its angular dependence in 
momentum space (Kohn singularity) results in the long-range 
Friedel oscillations of the two-particle potential in real space, 
inverting the sigh of the effective interaction. A diagrammatic 
expansion of the effective vertex at low fermionic density is 
known to result in the effective attraction in the second order 
of perturbation theory in 3D (third order in 2D) fsC 
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We apply this idea to the mixture of two types of spin- 
polarized cold gases with different masses and densities. Such 
two-component fermion systems have been realized either as 
mixtures of different atomic states of ^ Li ll37ll and '^^ K 113 811 
and also mixtures of ^Li-'^'^K ll39ll . For short-range interaction 
potential, the antisymmetry of the fermionic wave function 
results in the absence of diagonal intra-component terms in 




FIG. I. Effective interaction between fermions of the type A medi- 
ated by fermions of the type B in second order perturbation theory. 



the interaction Hamiltonian. The leading repulsive bare inter- 
action is reduced to the isotropic s-wave scattering between 
atoms of different kind with the scattering amplitude aq > 0. 
Such a system can be described by the following Hamiltonian: 

i=A,B p 

pp'q 

where i and j indicate species of two different types (A and 
B). The effective interaction Taa between fermions of the 
type A in the second order of perturbation theory (the low- 
est order where the effect exists) is mediated by the fermions 
of the type B and is represented diagrammatically in the 
Fig. [T] In the low-density limit ropp ^ 1 where ro 
is the range of the potential, the effective interaction term 
Fig. □ becomes T AAiqA) = -X'^nBsiqA), where qa 
Pfa-\/2(1 — cos 9) is the difference between the incoming 
and outgoing momenta and k/i with 9 the angle between 
them, A = 2ao y^p faP fb / is the modified 3D gas param- 
eter in the case of inter-component scattering. Furthermore, 

UBsiqA) - J {pB + ^)G {pB - ^) is the po- 

larization operator that corresponds to the Green's function 
loop in the Fig. [T] Integrating TAA{qA) at the Fermi sur- 
face with eigenfunctions of the angular momentum (Legendre 
polynomials Pi {cos{9)) in 3D), one can see that the most neg- 
ative eigenvalue f^^^ = /f^^(g^(6»))F,(cos(6'))^^ cor- 



responds to Z = 1 and the triplet superfluid state is formed in 
the system A. The critical temperature is then given by; 
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where Nfi and e Fi are the density of states at the Fermi sur- 
face and the Fermi energy of the species i. 

The static polirization operator in 3D is given by: 
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FIG. 2. Triplet critical temperature of the majority species A as a 
function of population imbalance 5 = pfa/pfb, at pfa > Pfb 
and fo ~ 0.8 in 2D. The minority component B remains in the nor- 
mal state in the second order of perturbation theory. Inset: The criti- 
cal temperatures for the species A and B in 3D at A ~ 1.9. 



then the critical temperature T^.^ is 
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is a dimensionless function of the population imbalance S = 
Pfa/pfb L30i,i4ft1 . One can see that f{S), and consequently 
T^^^iS) increases at small imbalance S >l, has a maximum, 
and then drops to zero at S ^ 1. Two factors contribute to 
the effect: on one hand, Kohn singularity increases with S and 
this enhances the effective interaction of the fermions of the 
type A; on the other hand, the density of the species B de- 
creases as 6 increases, and this eventually results in the total 
absence of pairing at S oo. The effective interaction for 
the species B is, in turn, mediated by the subsystem A. To 
obtain the critical p-wave temperature Tj:]^ one should swap 
indices A and B in the Eqs. (O - (|5]l. Then, the temperature 
T^'^j' decreases monotonically as d increases. At S = 1 the 
system is equivalent to the non-polarized 3D Fermi gas, and 
the superfluidity occurs in both subsystems simultaneously at 



Thus, at a finite imbalance 

(1) / / rrW 



n ' ^ £;^exp[-13/A2] il 

S > I there is a temperature window T^^ < T < T^^ where 
only the component A is superfluid, while the other compo- 
nent B remains in the normal state (See Fig. |2]i. 

Similar arguments apply to the 2D case, where cos lip, the 
2D eigenfunctions of the angular momentum, replace the Leg- 
endre polynomials, and qa = pfa^/'2{1 — cos(/3). The polar- 
ization operator 
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explicitly depends on ip (and thus on qa) only if pf A > Pfb- 
In the opposite case, ppA < Pfb, the angular dependence is 



completely excluded from the real part of the Eq. (O and all 
the eigenvalues Taa vanish. The population imbalance be- 
comes the necessary condition for the pairing to exist in the 
second order of perturbation theory. The p-wave superfluid 
state is then formed by the species with larger value of the 
Fermi momentum (subsystem A), with the critical tempera- 
ture given by the equation i30ll : 



(1) 



r^^^exp HV(4(<5 - l)/o')] 
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where fo is the 2D gas parameter that, as usual, is obtained 
by the renormalization of the two-particle vertex and that in 
2D contains a logarithm fo sa —l/2\nroPF- Just like the 
3D case, the critical temperature (|7]i has a maximum at some 
S 2 (See Fig. |2]l. The other, low-density subsystem B 
remains in the normal state at all temperatures in this order. 
Note that, from Eq. Q, T^^^ = at (5 = 1 in the second 
order of perturbation theory, and one needs to include the next, 
third-order terms of the expansion. To this order, the critical 
temperatures for both subsystems (note that for S ~ 1 the 
components A and B should behave similarly) are given by 
(4) Ti^) ^ ep cxp [-1/(6.1/3)] |31j]. 

Note that, according to Eq. (|2| and Fig. |2] at the gas pa- 



rameter fo ~ 0.8, the 2D maximum T 



(1) 



0.2eF is well in- 



side the currently experimentally achievable range (4)}]. The 
3D critical temperature is within the experimental range at 
A ^ 1.9 in Fig. |2] The applicability of the perturbation 
theory at these values of the gas parameter requires further 
justification, in particular, higher order terms in the perturba- 
tion theory which we leave for future work. Here we, how- 
ever, remind that, as was demonstrated in Ref. (3d], the zero- 
imbalance critical temperature at the gas parameter A ~ 1.3 
(same as that in the superfluid ^He) has the right order of mag- 
nitude ri^^ - IQ-^eF- 

In the present proposal we ignore possible off-diagonal ex- 
change terms in the Hamiltonian ([T]) that scatter Cooper pairs 
between different bands in the two-band model 142(1 . These 
terms will create a superfluid state in both bands simultane- 
ously 1I40I1 . In the present case where components A and B 
represent two types of atoms (e.g., ^Li and ^"^K) such off- 
diagonal terms cannot exist. In addition, we see that when 
both species are completely polarized, the critical tempera- 
tures splitting is mostly affected by the population imbalance, 
or the ratio of Fermi momenta of different components, but 
not the by the ratio of their masses. 

RF spectroscopy based detection of Majorana fermions: 
Topological superfluids of the type considered here are char- 
acterized in Z? = 2 by Majorana fermion excitations that are 
localized in defects of the superfluid order parameter. Such 
zero-energy states can be detected as resonances in tunnel- 
ing near zero energy, which can be realized in cold atomic 
system using RF spectroscopy Let us consider 

a two-component mixture of spin-polarized ^'^K as the ma- 
jority species A and ''Li atoms as the minority species B. 
We first note that for the purpose of imaging MFs using RF 
spectroscopy, the fermions in the components A and B can 
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be discussed independently because of the absence of inter- 
conversion between the two species of atoms. As discussed 
above, from the Kohn-Luttinger effect, the majority A atoms 
can be expected to be in a topological superfluid phase, while, 
in D — 2 the minority component B remains uncondensed. 
The uncondensed B atoms, however, pose no problem for the 
detection of the MFs in component A as we clarify later. Sup- 
pose the ^"^K atoms in the superfluid are in the 4^51/2 hyper- 
fine ground state, \i) = |F = 9/2, mp = -7/2). In this 
case, the application of a two-photon Raman pulse can trans- 
fer the '^"K atoms from the state |i) to another hyperfine state 
\j) ^ \F ^ 1/2, mp = -5/2). Within the rotating wave- 
approximation (RWA), the frequency difference between the 
two Raman lasers Vl can be considered to be the detuning be- 
tween the energy levels of the internal states \i) and The 
effective Hamiltonian describing such a pair of Raman beams 
within the RWA is written as 

HRan^an^ I dr {/(r ) (r ) V, (r) + /i.c] + 1] (t") V', )} , 



(8) 

where ij^l ^ (r) are the creation operators for atoms in states |i) 
and \ j) respectively. Here we assume that a localized Raman 
beam with a beam waist width w « 1.5 /xm is much smaller 
than the typical distance 10/im between vortices [,46.1 is being 
used. The function f{r) is localized within the beam waist. 

The above Hamiltonian Hfiaman transfers atoms from the 
superfluid of ""'K atoms in the \i) state to atoms in the \i) 
state, which can be thought of as being in an empty reservoir. 
This is similar to a tunneling probe in a superconductor, where 
electrons can tunnel from the superconductor into the leads. 
For topological superfluids without vortices and other defects, 
which have a unique ground state 1 4*) together with a finite 
quasiparticle gap Aq to the lowest excited state l^*!), the rate 
of generation of atoms in state |j) must vanish for energies 

< Ep — Ao (where Ep is the fermi-energy of the super- 
fluid of \i) atoms) because of energy conservation. In contrast, 
topological superfluids with vortices contain zero-energy MP 
excitations bound to vortices. The MFs are associated with 
self-hermitean creation fermionic operators 7, which com- 
mutes with the total mean-field Hamiltonian for the topologi- 
cal superfluid. As a result if |^) is a lowest-energy state of the 
Hamiltonian, = 7!^') is another state with exactly the 
same energy (measured relative to the fermi-energy). There- 
fore the ground-state is degenerate in this case and it is possi- 
ble to generate particles in the state \ j) even when the Raman 
beams are tuned to be on resonance at = Ep . In fact, the 
rate of generation of fermions in the state \ j) will have a peak 
at n = Ep given by 



R = Y,\{'i'n\HHa 



^■^)\^s{E - E„ - n-Ep)(x s{n) 

(9) 

where I'ifn) are excited states with energy E„ and E is the en- 
ergy of the ground state j^I^). Note that, while the absorption 
of a Raman beam transfers the many -body state from 1 4*) to 
IvPi) in principle, this does not lead to any observable time- 



dependence of the particle-generation rate R since the matrix 
elements Hjiaman for localized Raman beams are identical 
for the states |4') and |^i). Therefore, the rate of particle- 
generation R is independent of which initial topologically de- 
generate state I'if) or I'ifi) one starts with, consistent with the 
idea of local non-observability of the topological degree of 
freedom. The tunneling operation induced by the Raman- 
pulse cannot be used to read-out the MP degree of freedom 
unless the two vortices are brought together to induce a split- 
ting - effectively fusing the vortices. 

Another interesting feature of using a mixture of fermions 
containing two different species of fermions is the topolog- 
ical robustness despite the presence of an even number of 
fermi surfaces. This is typically a problem in spinful chiral 
p-wave superconductors such as SrRu04 where a vortex con- 
tains a spinful Majorana fermion with a spin-degeneracy. In 
that case, transverse magnetic fields which mix the two spin- 
species would split the Majorana fermions and lift the topo- 
logical degeneracy. The observation of true MFs in this sys- 
tem requires the use of half-quantum vortices l47ll48ll . which 
are effectively vortices in one of the spin species to trap Ma- 
jorana fermions without a spin degeneracy. The cold atomic 
system proposed here provides a unique robustness in this re- 
gard, because here the two species or pseudo-spins are really 
two different elements Li and K in this case, so that the Majo- 
rana fermions are protected by the additional number conser- 
vation symmetry in this system. Therefore, even if one used a 
regular vortex, which produced a phase winding in each of the 
species separately, there would be a pair of Majorana fermions 
in each vortex, which cannot be split simply because it is not 
possible to convert an atom of one element to another A simi- 
lar protection would continue to exist if one of the species (say 
Li) was not condensed, as is the case in 2D (Pig. 2). Conven- 
tionally, one could worry that the thermal Li gas serves as a 
fermion reservoir which is known to destroy the fermion par- 
ity protection of Majorana fermions. However, since Li atoms 
cannot convert into K atoms, the Majorana fermions associ- 
ated with the K topological superfluid remain robust. 

In D = 3, the topological superfluids considered here are 
characterized by topologically protected gapless Dirac points 
in the bulk (Weyl fermions) and flat Majorana modes on the 
surface y5|-|7|]. Such momentum resolved gapless spectrum 
can in principle be detected by recently developed momentum 
resolved photo-emission spectroscopy for cold atoms [41, 49)]. 

Conclusion: Using Kohn-Luttinger effect we show that in a 
mixture of two species of individually spin-polarized fermions 
with population imbalance it is possible to create topologi- 
cally non-trivial superfluid states characterized by Majorana 
fermion excitations at order parameter defects. The Kohn- 
Luttinger critical temperatures of the individual superfluids 
are split as a function of the population imbalance and are 
well inside experimentally achievable range in current exper- 
iments. The Majorana fermions in the individual components 
of the superfluid can be detected by RP spectroscopy. Be- 
cause of the absence of inter-conversion of the different types 
of atoms, the Majorana fermions are detectable even if one 
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of the components remain uncondensed. Our proposal to cre- 
ate topological spin- triplet ^j- wave superfluids using the Kohn- 
Luttinger effect, coupled with the methods to observe the non- 
trivial topological excitations, gives a complete description of 
a promising new way to create and analyze a topologically 
non-trivial superfluid in both two and three dimensions. 
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